Integer quantum Hall effect of interacting electrons: dynamical scaling and critical 

conductivity 
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We report on a study of interaction effects on the polarization of a disordered two-dimensional 
electron system in a strong magnetic field. Treating the Coulomb interaction within the time- 
dependent Hartree-Fock approximation we find numerical evidence for dynamical scaling with a 
dynamical critical exponent z = 1 at the integer quantum Hall plateau transition in the lowest 
Landau level. Within the numerical accuracy of our data the conductivity at the transition and the 
anomalous diffusion exponent are given by the values for non-interacting electrons, independent of 
the strength of the interaction. 
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While the occurrence of quantized plateaus in the Hall 
conductivity of two-dimensional systems in strong mag- 
netic fields is the most striking aspect of the quantum 
Hall effect (QHE) §, it is the transitions between these 
plateaus that have recently attracted a lot of attention, 
both experimentally and theoretically Q. Our under- 
standing of the QHE is not complete without an under- 
standing of the plateau transitions. Furthermore, the 
plateau transitions are among the most extensively stud- 
ied examples of quantum phase transitions. The origin of 
the plateaus is the localization of the electrons due to dis- 
order, with the plateau transitions corresponding to the 
localization-delocalization transitions in the disorder po- 
tential Q . Theoretical studies have focused on the effect 
of disorder on non-interacting electrons in a strong mag- 
netic field. From numerical and analytical calculations 
the following picture emerged for weak disorder Q : in ev- 
ery Landau level exists a critical energy E c at which the 
localization length £ diverges as a power law \E — E c \~ v 
with a localization length exponent v — 2.35 ± 0.03 in- 
dependent of Landau level index and correlation length 
of the disorder potential ||. At all other energies elec- 
trons are exponentially localized. In the region of local- 
ized states the longitudinal conductivity vanishes and the 
Hall conductivity is quantized in multiples of e 2 /h. At 
the critical energies the Hall conductivity jumps by e 2 jh 
and the longitudinal conductivity takes on a finite value 
a c with <j c w 0.5e 2 //i in the lowest Landau level 
The eigenstates at the critical energies show multifractal 
fluctuations leading to correlation functions being char- 
acterized by an anomalous diffusion exponent r/ w 0.4 [Q. 
Dynamical correlation functions show scaling behavior as 
a function of the variable qL u , where q is the wave- vector 
and L u = 1/ V pqHlu is the relevant frequency-dependent 
length scale, the effective system size with mean level 
spacing Hu> Q|. Comparing this result to the definition 
of the dynamical critical exponent L w oc uj^ 1 ^ and ob- 



serving that the density of states (DOS) po is finite at 
the transition, the dynamical critical exponent is found 
to be z = 2. 

Experimentally, strong evidence for scaling behavior 
near the plateau transitions was observed Q. In par- 
ticular, the observed values for the localization length 
exponent v — 2.3 ± 0.1 [|| and the critical conductiv- 
ity <j c w 0.5e 2 /h are in remarkable agreement with 
the theory for non-interacting electrons ||. However, 
in high-frequency measurements scaling with a dynami- 
cal critical exponent z = 1 was observed in contrast to 
the value z — 2 for non-interacting electrons ||. This 
discrepancy has been attributed to the influence of the 
Coulomb interactions between the electrons that are al- 
ways present in the experiments JlO[ ] . In fact, numeri- 
cal simulations in which the interactions are treated in 
self-consistent Hartree-Fock (HF) approximation show a 
linear suppression of the DOS at the Fermi energy irre- 
spective of the position of the Fermi energy jll],[l2| ■ Ac- 
cording to the argument presented above at the critical 
energy this non-critical effect leads to a reduction of the 
dynamical critical exponent z from 2 to 1. The influence 
of the Coulomb interactions on other critical properties 
of the QH transitions was also studied within the self- 
consistent HF approximation [fl2|| . From the scaling of 
the participation ratio of the HF eigenstates it was found 
that the localization length exponent v and the fractal di- 
mension D2 = 2 — 7/ were not changed by the interactions. 
In lieu of the critical conductivity the Thouless numbers 
of the HF eigenenergies were studied and also found to 
be unchanged by the interactions. However, not much is 
known about the relation of the conductivity of an inter- 
acting electron system and the Thouless numbers of the 
HF eigenvalues. 

In this paper we overcome some of the deficiencies 
of the previous studies by directly calculating the irre- 
ducible dynamical polarization H ln '(q,uj). This allows us 
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to obtain the critical quantities z, a c , and r\ from the 
scaling limit of the polarization. Let us first sketch our 
approach. As the starting point for our numerical cal- 
culations we use the self-consistent HF approximation 
JTT| , fL2"| . The corresponding conserving approximation for 
the two-particle Green's function is the time-dependent 
Hartree-Fock (TDHF) approximation Q. Within this 
approximation we evaluated the dynamical polarization 
II lrr (g,a;) including vertex corrections. The disorder- 
averaged polarization is extrapolated to infinite system 
size and vanishing q and lj. In this limit, the polarization 
obeys the scaling form 



n irr ( 9 ,c) 



a * q 2/ e 2 x m 



IUJ 



(1) 



The Onsager conductivity a* is related by the Einstein 
relation 





0.7 r 




0.6 




0.5 


o 




Q. 


0.4 


ES 


0.3 






0.2 




0.1 












0.1 

UyL 



100 
169 

256 
400 
625 
900 



0.2 



FIG. 1. Single-particle density of states at the critical en- 
ergy as a function of for N<s> between 100 and 900 and 7 
between 0.15 and 1. 
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to the diffusion coefficient D(q,uj) and the static suscep- 
tibility Xq 1 ■ F° r vanishing wave- vector it coincides with 
the Kubo conductivity a — e 2 / iq 2 . Due to the oc- 
currence of the Coulomb gap even at criticality the sus- 
ceptibility x 1 " cx q. The conductivity a*(q, ui) is found to 
be a function of the scaling variable x = q 2 /x 1 "^ with 
the limiting behavior 



a*(x) 
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Since x oc q/u> the Coulomb gap changes the dynamical 
critical exponent z from its non- interacting value 2 to 1. 
On the other hand, we find that the critical DC conduc- 
tivity a* = 0.5 ± 0.1e 2 /h and the anomalous diffusion 
exponent 7/ = 0.4 ±0.1 are independent of the strength 
of the interactions and in agreement with previous calcu- 
lations for non-interacting electrons (J,|^,|lJ]. These are 
the central results of our study. 

In the remainder of the paper we present our model, 
our calculations and our results, relegating a more de- 
tailed discussion to a separate publication (T5j . We con- 
sider the situation where the cyclotron energy fiuj c is 
much larger than both the disorder and the Coulomb en- 
ergy and the Fermi energy is in the lowest Landau level 
onto which we project the Hamiltonian. As a model for 
the disorder we use the random Landau matrix fl6f| . The 
strength of the disorder potential is characterized by the 
level broadening in self-consistent Born approximation T 
and its range is taken to be zero. The electrons inter- 
act via the bare Coulomb interaction V(q) = 2ne 2 / nq. 
The strength of the interactions is given by the ratio 
7 = (e 2 /nl)/T, where k is the dielectric constant and 
I is the magnetic length. The square systems of linear 
size L with periodic boundary conditions contained N 
electrons and from N& = L 2 /2irl 2 = 169 to 900 flux 



quanta. In the presence of interactions physical quanti- 
ties depend on the filling factor vp = N/N&. We choose 
vp = 1/2 in all our calculations so that the Fermi energy 
coincides with the critical energy of the plateau transi- 
tion. Including the HF interaction the eigenvalues e a and 
eigenf unctio ns \a) of the Hamiltonian are calculated self- 
consistently | Q , jl2fl . A major simplification arises in this 
calculation since in the lowest Landau level the exchange 
interaction 

V cx (q, q') = -Sq^q,V2^e 2 l K - 1 e^ l2 /%(q 2 l 2 /4) (4) 

and hence the HF interaction 

V UF (q, q') = & q ,-f (V(q) + V cx (q)) (5) 

become local (Io is the Bessel function) Jlj],|l8| ■ 

Perhaps the most characteristic feature of the HF spec- 
trum is the formation of a Coulomb gap in the single- 
particle DOS p9| . In particular, the DOS at the Fermi 
energy p{Ep) decreases linearly with 1/jL. This is shown 
in Fig. The suppression of the DOS is not re- 

lated to the phase transition but occurs for all values 
of the Fermi energy. Only when the Fermi energy coin- 
cides with the critical energy does it change the critical 
dynamics of the system. 

The bare polarization n°(q, q', us) is given by 

n°(q,q>)=^(a|/?) q (/3|a) q , 

<x,P 

U±-ff}) Ml -/a) 



huj + e/3 — e a — iS fiuj + 6/3 — e a + id 



(6) 



where (/3|a) q = (/3|e _lqr |a), f a is the Fermi function 
(we take T = in our calculations) and 6 is an in- 
finitesimal parameter of the order of the mean level 
spacing. After averaging over disorder (we use between 
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FIG. 2. Feynman diagrams for the irreducible vertex (top) 
and the screened polarization (bottom) in time-dependent 
Hartree-Fock approximation. 



13 and 100 realizations of the disorder) II becomes 
translationally and rotationally invariant II (q, q',w) — 
<5(q + q')II (g, uS). The real part of the static bare po- 
larization coincides in the limit q — ► with the single 
particle DOS, lim^o ReIl (q,LO = 0) = Xa = p( e f) |j§. 

The TDHF approximation for the irreducible and the 
screened polarizations are shown in Fig. ^. In the lowest 
Landau level the ladder sum for the vertex corrections 
can be transformed into a RPA-like bubble sum with the 
Coulomb potential replaced by V ex (q) pc| , [l5|| . Similarly, 
the screened polarization II scr is related to the bare po- 
larization through a bubble sum containing the HF po- 
tential Vhf(c?) @, 

n h ' r (q, U ) = n°(q, w) (1 - VUq)n irr (q co)) , (7a) 
n scr (q,w) = n°(q,^) (1 - V HF (q)n SCI '(q,w)) . (7b) 

In principle, the vertex corrections should be applied to 
IT (q, q') since it becomes diagonal only after disorder 
averaging. We checked for small system sizes that using 
Eq. (Q) only introduces errors comparable to the statis- 
tical uncertainties for the relevant parameter range. 

Since the approximations leading to II lrr and IP cr are 
conserving and due to Eq. (Q), the polarizations can be 
parameterized in the limit of infinite system size and van- 
ishing q and co through a susceptibility Xq and a diffusion 
coefficient D(q,Lo), 



Tl A (q,Lo) 



V 



D A q 2 



(8) 



with A — {0, irr, scr}. The Einstein relation a (q,u) = 
e 2 XqD A (q 1 lo) defines the Onsager conductivity 0* . From 
Eq. (0) it is apparent that the dressing of Xq an d D leaves 
the Onsager conductivity unchanged 



e 2 X°D° 



e 2 Xq TT D iTr 



,scr^scr 



(9) 



q-^ " /Vq & Xq 

For vanishing wave- vector the Onsager conductivity coin- 
cides with the Kubo conductivity cr*(q — 0,lu) — <j{lo) = 
lim^o e 2 LoU 1T1 '(q, lo) /iq 2 - 

The polarization can be rewritten in terms of a* , 



W"(q,Lu) = x q 



ax 



ie 2 /%' 



(10) 



with x = q 2 /Xq T fcw- If dynamical scaling holds then a*x 
is a function of q z /lo. In order to extract the scaling be- 
havior from our numerical data for finite system sizes, 
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FIG. 3. Static irreducible susceptibility for N<s = 900 and 
7 = 0.2, 0.4 and 0.6. The finite intercept at q — scales like 
l/L to zero. 
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FIG. 4. Extrapolation of ImII lrr (q, tj)/xlj rr as described in 
the text. x/2ir = 0.5, 1, 1.5, 2, 4, 6 (bottom to top). 



we first calculate the static susceptibility x 1 " from the 
real part of the irreducible polarization. Due to the oc- 
currence of the linear Coulomb gap (cf. Fig. , x 1 " °^ 1 
for small q in the thermodynamic limit, as is shown in 
Fig. & 

Next, the imaginary part of IT" is divided by the nu- 
merically obtained Xq T - The resulting function is then 
simultaneously extrapolated to infinite system size and q 
and lo to zero for fixed value of x as is shown in Fig. |] 
|^,^5|. This choice of scaling variable allows for reason- 
able extrapolations, while other choices like q 2 / pohco do 
not give satisfactory results. Since x irr is asymptotically 
linear in q, x oc q/co in the thermodynamic limit and the 
dynamical critical exponent z = 1. 

Figure [| shows the resulting conductivity a* as a func- 
tion of x. In the limit x — > 0, a* coincides with the DC 
Kubo conductivity a c . We find that a c does not depend 
on the strength of the interaction up to 7 = 0.6 and is 



3 



2.0 
1.5 

1.0 

0.7 

05, 

0.7 

0.5, 

Y 

0.7 
0.5 




1 10 

xll% 

FIG. 5. Conductivity a* as a function of the scaling vari- 
able x = q 2 /Xq T Tiuj for 7 = 0.2, 0.4 and 0.6 (bottom to top). 
The curves are shifted vertically. 



given by a c — 0.5 ± e 2 jh. In the opposite limit x — * 00, 
the conductivity reflects multifractal density fluctuations 
and scales like x 1 ^ n A numerically more accurate 

determination of r] is possible from the scaling of the 
static polarization |15| |. We find 77 = 0.4 ±0.1 indepen- 
dent of the interactions and in agreement with published 
data for non- interacting electrons p]|l4|]. Comparing our 
result to |Q we see that the irreducible polarization re- 
flects the multifractal exponent D2 of the HF cigenstates 
similar to the case of non-interacting electrons. 

Let us make a few remarks on our results. Although 
the conductivity a c is finite at the transition, the van- 
ishing susceptibility x 1 " implies non-diffusive behavior 
of the HF-particles, D(q,u) oc q^ 1 for q — > 0. A cen- 
tral issue is the validity of the HF approach that has 
recently been questioned |pl[| . The change in the dynam- 
ical critical exponent is due to the use of the unscreened 
Coulomb interaction. The dielectric function k is given 
by K(q,w) = 1 + ^(<j)n irr (g,w). From eq. (§) it follows 
that screening behaves differently in the limits of small 
and large x. In the limit of small x, that is relevant to 
transport, screening is inefficient and k — ► 1. The use of 
the unscreened HF approximation thus might be justified 
for the calculation of the frequency-dependent conductiv- 
ity. On the other hand, in the limit of large x we observe 
static screening, k{c[,lj) = 1 + Xq r V(<l)- For slow pro- 
cesses, during which the ground state has time to relax, 
the statically screened HF approximation is more appro- 
priate. It leads to a finite susceptibility and is described 
by z = 2 |lq ]. A final judgment on the validity of the HF 
approach to the critical quantum Hall system might only 
by possible if one succeeds in a renormalization group 
treatment of the neglected correlations |h],^2| . 

In summary, we have presented results on the influ- 
ence of Coulomb interactions on the critical properties 
of the integer quantum Hall system. Within the time- 



dependent Hartree-Fock approximation we find the criti- 
cal conductivity a c and the anomalous diffusion exponent 
r\ are unchanged by interactions while the dynamical crit- 
ical exponent z changes from 2 to 1. The origin of this be- 
havior is the non-critical reduction of the single-particle 
density of states. 

This work was performed within the research program 
of the SFB 341 of the DFG. 
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